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Abstract: The application of a high purity Germanium (HPGe) 7 spectrometer in determining the fuel el- 
ement burnup in a future reactor is studied. The HPGe detector is exposed by a ^^Co source with varying 
irradiation rate from 10 xlO'^s"^ to 150 xlO'^s"^ to simulate the input counting rate in real reactor envi- 
^ ronment. A ^''^Cs and a ^^^Eu source are positioned at given distances to generate certain event rate in the 

■ detector with the former being proposed as a labeling nuclide to measure the burnup of fuel element. It is 

shown that both the energy resolution slightly increasing with the irradiation rate and the passthrough rate at 
high irradiation level match the requirement of the real application. The influence of the background is studied 
in different parameter sets used in the particularly developed procedure of the background subtraction. It is 
C/3 ■ demonstrated that with the typical input irradiation rate and '^^'^Cs intensity relevant to deep burnup situa- 

tion, the precision of the ^^'^Cs counting rate in the current experiment is consistently below 2.8%, indicating 
^ ■ a promising feasibility of utilizing an HPGe detector in the burnup measurement in future bed-like reactor. 
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^ . 1 Introduction 

o 

' The fuel balls undergo a multipass cycle in the modular pebble bed reactors (MPBR). Because of the serious 
• : inaccuracy of the computational method, which is conventional for performing the in-core fuel management in 
' the existing water reactor, a non-destructive determination of the burnup is desired to provide the distributed 
ff^ • controlling^stem with an online circulation/discharge judgement on a pebble- by-pebble basis in some bed- like 
, reactors[lJ,|3|. Gamma ray spectrometry has been proposed as an effective non-destructive method to determine 
' the burnup [3|-[7| as well as the spatial irradiation distribution [8, 9] by measuring the activity (activity ratio) 
of given monitor nuclide(s). Among all kinds of fission products, ^^''Cs has been proposed to be one of the 
^ , effective burnup indicators for its long life-time and the clear correspondence to the burnup of the fuel pebble 
with rather resistance to the power history 10, 11|. Nevertheless, the fission of the fissile materials in MBPR 
produces various radioactive nuclides that emit amount of 7 peaks over the whole spectrum, and hence, as 
shown by various simulation results [H, E2-14| and some testing results, the ^^^Cs full energy peak is likely 



obscured by the neighboring 7 rays in a typically short cooling time. For instance, in almost all the stage 
of the burnup, the peak of 658 keV from ^^Nb is presented. Thus, in order to identify the 7 peak from the 
^^^Cs 662 keV, a 7 spectrometer with energy resolution of better than 2 kcV is highly desired. In this regard, 
the High Purification Germanium (HPGe) 7 detector has been widely proposed in the determination of fuel 
pebble burnup as a non-destructive method. 

Another feature of the fuel element with a certain burnup in the certain future reactor is the high 7 ra- 
dioactivity. According to the prediction, the radiation exposing to the HPGe detector varies from tens of 
xlO^s"^ to more than 150 xlO^s"^ , depending on the shielding condition, the geometric collimation and the 
cooling time. In the design of the spectrometer, although the dead time of the spectrometer can be lowered 
by decreasing the irradiation level through a strict geometrical collimator, the net count rate of the labeling 
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^^^Cs is also suppressed and a larger statistical error is expected. Therefore for the MPBR under design, it is of 
extreme importance to optimize the working parameters, including the shaping time, geometrical collimation 
and shielding, to have the activity determined with a minimized experimental error limit within the typical 
real time (tens of seconds) after very short cooling time (typically tens of hours). 

In this paper, we present the main results in experimental conditioning of an HPGe 7 spectrometer. The 
passthrough curves at different shaping time and flattop are measured. For the appropriately chosen shaping 
time, the peak resolution is measured with various radiation background levels. Furthermore, the precision 
of measuring the net counts of a ^^'^Cs source is studied in detail by varying the real time, the net rate of 
the nuclide of interest (NOI) and the radiation background level. In Section II, the experimental setup and 
conditions are described. Section III presents the main results and Section IV is the summary. 

2 The experimental setup 

The whole spectrometer consists of an HPGe 7 detector, an electrical cooling device (X-Cooler II) and a 
suitable fast electronics (DSPEC+) from ORTEC. The cylindric HPGe crystal is 43 mm in height and 62.6 mm 
in diameter, with a (plO.Smm x 30mm hollow copper electrical pole fed in the axial center. The efficiency of 
the HPGe detector for ®°Co 7 rays is about 30% at the working high voltage +2300V. The preamplifier adopts 
the light feed-back technique to reduce the pile-up effect at high incident rate. 

Three 7 sources are used in the experiment. A ®°Co is applied to simulate the irradiation background from 
10 to 150 xlO^s"^ (dubbed from COl to C015 in the text) by varying its distance to the detector. The intense 
incident counting rate is controlled by the ^°Co source throughout the test unless specially specified. A ^^""Cs is 
placed at two positions in the vicinity of the detector to generate about 30 (CSl) or 120 (CS2) counting rate 
in the detector to mimic the counting rate of a low and a high burnup element, respectively. Correspondingly 
a ^^^Eu source is also introduced to investigate the performance of the spectrometer varying with the 7 energy. 
At each geometrical configuration, the measurement is done over 2 settings of the shaping time and the flattop 
(SHPl and SHP2) of the main amplifier for 1 or 10 groups, with 5 files in each group. Totally 2904 files are 
acquired. With these measurements it is possible to study the repeatability and the statistic fluctuation of 
multi observables. Table [T] summarizes the experimental conditions of the multi source experiment. 



Table 1. Conditions in the multi source experiment 





Shaping time 


^^''Cs Intensity rics (s ^ ) 


''"Co Irradiation 




and flattop Tj (/is) 


and real time (s) 


rate nco (xlO^s"^ ) 


Subscript 


SHPl SHP2 


CSl CS2 


COl ... C015 


Parameter 


T, = 0.8 T, = 1.2 


ncs=30 ncs = 130 


nco = 10 ... 7ico = 150 


Setting 


Tf = 0.8 rj = 0.6 


Tr = 10 = 25 





3 Results and discussions 

3.1 Passthrough and energy resolution 

Figure [T] presents the passthrough curves measured with SHPl and SHP2. The abscissa depicts the incident 
7 counting rate characterizing the radiation background level exposed on the detector. The low threshold is 
0.3 keV. It is shown that at low radiation level, the passthrough rate increases linearly with the input rate. At 
about 100 xlO^s"^ , however, the passthrough starts to deviate with different shaping time and undergoes a 
plateau up to an incident rate of about 200 xlO^s"^ . Above 200 xlO^s"^ , the passthrough starts to decrease 
with the input counting rate due to the rapidly increasing dead time. Unlike the situation with low input 
rate, where the passthrough shows insigniflcant dependence on the shaping time, the maximum passthrough 
on the plateau decreases with increasing the shaping time and flattop as shown in the flgure. Quantitatively 
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the passthrough is 10 xlO^s"^ higher with SHPl than that with SHP2. In the real application, the amplifier 
parameters and the geometrical coUimation are so designed that the system is running on the passthrough 
plateau to acquire the highest statistics while keeping adequate energy resolution since the 7 energy peaks in 
the proximity of the ^^^Cs indicator enrich. As shown later, one can not switch the shaping time down to very 
small value for the rapidly degraded energy resolution. 




Fig. 1. (Color online) Passthrough rate as a function of incident counting rate with SHPl and SHP2. The 
solid, dashed and dotted lines represent certain dead time of 0%, 30% and 60%, respectively. 



Figure [5] depicts the dependence of FWHM on the shaping time for ^^^Cs and ""Co peaks. With increasing 
the shaping time, the FWHM of the full energy peaks decreases rapidity below about 4 fis and then saturates 
at 1.3 keV and 1.8 keV for ^''^Cs and ''''Co full energy peaks, respectively. Shaping time of larger than 0.8 fiS 
is preferred with respect to the required resolution of 1.8 keV for ^^^Cs. But because of the high expectation 
of passthrough, optimization is required. In this paper, only the results with SHPl and SHP2 are relevantly 
presented since these two conditions meet the requirements of the design of 7 spectrometry for the future real 
application. 
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Fig. 2. (Color online) FWHM of the full energy peaks as a function of shaping time for the ^"^^Cs peak (662 
keV) and the '^'^Co peaks (1173 and 1332 keV). 



The energy dependence of FWHM is plotted in Figure [3] at two typical radiation background levels 30 and 
100 xlO^s"^ . The data are fitted with a formulae of ay/E^+b as shown by the dashed curves. Although on 
average the data points follow the curves well at both irradiation rates, rather fluctuation of the data points 
is visible. This is attributed to the relative intensities at different energies of the ^^^Eu peaks. The full energy 
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peaks with lower intensities, for instance, at =411 or 1090 keV, suffer more from the background level 
induced by the Compton plateau of the intense ®°Co source. And equivalently, the FWHMs become smaller if 
the incident rate is lowered to 30 xlO^s"^ as the Compton background decreases. As expected, the FWHMs 
of the two ^"Co peaks at 1173 and 1332 keV are not changed by the radiation rate. 
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Fig. 3. (Color online) FWHM of the full energy peaks as a function of peak energy at different input rate. 



FigureSlpresents the FWHM of the ^^^Cs full energy peak as a function of incident counting rate at different 
shaping time configuration SHPl and SHP2, respectively. Specially the irradiation rate corresponding to the 
left panel is controlled by moving the ^^^Cs itself. In this case the FWHM keeps constant with the input 
counting rate. However, in the test where the input counting rate is controlled by moving the ""Co while the 
^^'^Cs position is fixed, the FWHM exhibits an increasing trend with the input counting rate. At typical input 
rate of 30 xlO^s"^ corresponding to deep burnup, the FWHM of ^^^Cs full energy peak is below 1.8 keV. The 
energy resolution of ^^^Cs peak is degraded by less than 20% at 150 xlO^s^^ input rate. 
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Fig. 4. (Color online) FWHM of the full energy peak of ^'^^Cs as a function of input counting rate at different 
shaping time for the situations where the incident counting rate is controlled by ^"^^Cs itself (a) and by the 
'^''Co source (b). 



3.2 Precision of the net counts of the NOI 

The relationship between the activity of the source A and the net counting rate of the full energy peak n 
recorded by the detector is written as 

n = ATe (1) 
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where e is the full energy detection efBciency and F is the absolute strength of the peak, n is obtained by 
subtracting the background B from the total counting T in the range of interest (ROI) corresponding to the 
full energy peak range via 

n= — ^ 2 

t t ' ' 

where N is the net counts recorded in the live time t of measurement. The standard deviation of the net 
count cTjv writes 

On = ^/{otY + {(TbY ^VN + 2B (3) 

In real application, the total background B is derived by extrapolating the integrated counts B^ in two 
predefined windows (with CH,, channels for each) at the board to the whole ROI (with 2CHhr channels in 
total) via 

^-^-"^ w 

It is found in the analysis that the net counting rate depends slightly on the range of ROI for the same 
spectrum. In order to minimize the fluctuation caused by the change of ROI range, we fix the number of 
channels of the ROI for each given energy peak at a certain shaping time. Figure [5] shows the sketch of the 
procedure of the background subtraction. On both sides of the ROI, each of which is CH^r to the center of 
the peak, a fixed window with CHi, channels at the broad is selected to derive the background by a linear 
fitting to the data points in these two windows (the green shadowed area). The net peak (the cross symbols) 
is consequently obtained by subtracting the fitting background(the dashed line) from the total spectrum (the 
solid histogram). The net area of the peak is then obtained by fitting the net peak with a Gaussian function 
(the solid curve). 
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Fig. 5. (Color online) Procedure of the background subtraction. The histogram shows the original measured 
7 spectrum, and the green shadowing area depicts the windows used to fit the background with CHt channels 
for each. The distance between each broad of the ROI to the peak is denoted as CHhr- The cross symbols 
represent the net peak after subtracting the fitted background denoted with the dashed line, and the pink 
curve is a Gaussian fit to the net peak. 



Writing the averaged count of the background in a single channel as b for a given radiation background 
level, from Formulae ([3]), one then gets the uncertainty of net area 

/ CH^ I CH CH^ I 

'yM = jT + al^^ ^Jn + B^ + B^ = V ^ + bCH^l + CH,,JCH,) (5) 

It is shown that the uncertainty of the net area or net rate depends on the background level b as well as the 
net area N . Under a given radiation condition, the contribution of the background to ctjv decreases with the 
width of the background window because the fluctuation of the background level, shown with the dashed line 
in Figure O is weakened by averaging on more channels, as indicated by Formulae ([5]). Figure [S] presents the 



No. X 



6 



distribution of the background (the left panels) and the net area (the right panels) of ^^^Cs peak at different 
background window width CHi, for CSl at about 30 xlO^s"^ incident counting rate. With CHb tuning from 1 
to 9, the distributions of the net counts and the background become narrower, indicating a decreasing standard 
deviation with CHt,. 
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Fig. 6. (Color online) Distributions of the background B (left) and the net counts of ^^^Cs peak A'^ with 
different number of channels CHb of the background window (see Figure [5]) used to subtract the background 



Figure [7] presents quantitatively the relative standard deviation of the background <Jb/B (upper) and of the 
net area /N (lower) of ^^^Cs peak as a function of CHb for CSf and CS2 at 30, 50 and 100 x lO^s"^ incident 
counting rate, respectively. The curves are the fitting with formulae -^/ci + C2 /Nb with Ci and C2 being two 
fitting parameters. The relative standard deviation of the background is higher at lower irradiation rate while 
the one of the net area decreases with lowering the radiation background level. The effect of counting time is 
also clear by comparing the group of CSl and CS2. With longer measuring time (CS2), both the background 
and the net area exhibit less uncertainty. At the typical radiation level of 30 x lO^s^^ incident rate for the deep 
burnup fuel element with a ^^^Cs counting rate at about 130 s"^ , the precision of the net rate is better than 
3%. 

Figure m shows the relative standard deviation profile of the net area of ^^^Cs peak as a function of ®°Co ir- 
radiation rate for CSl and CS2, respectively. The symbols are calculated via Formulae ([3]) with the total area 
T and the background B. The curves , on the other hand, are derived from the corresponding root mean square 
(RMS) of the net area distributions, as shown exemplarily in the right columns of Figure[6]at the corresponding 
radiation level by RMS/N . It is shown that the curves follow the symbols well with few exceptions where the 
measurements are repeated for less than 5 times. It is shown that the uncertainty of N increases near linearly 
with the radiation background level, showing similar trends for CSl and CS2. For the group of CS2, for which 
the ^^^Cs source is closer to the detector and the real time is longer, the uncertainty is much smaller. 

More investigations including the behavior of ^^^Eu peaks in the same measurements reveal more information 
on the influence of the background level. Figure |9] shows the correlation between the relative standard deviation 
for all energy peaks below ®°Co 1173 keV as a function of the ratio of the net rate and the irradiation rate. 
Again the relative uncertainty is calculated via Formulae Q. The pink area corresponds to the group CSl and 
the red area to the group CS2. It is clear that by increasing the measure time and moving the source closer 
to the detector, the overall uncertainties are reduced. Since the background level induced by ^'^Co varies with 
the energy and the net rate differs among the ^^^Eu peaks, the uncertainty exhibits a rather broadening at a 



No. X 



7 




given scaled net rate n/ 1. The open triangles and the open squares with quadratic fitting (the dashed curves) 
present the evolution of the uncertainty of ^^''Cs peak for CSl and CS2, respectively. From this curve we can 
estimate the precision of measuring the ^^^Cs net rate if knowing roughly the ratio of the net counting rate 
and radiation background level. By neglecting the detailed distribution of the background sources, a precision 
of 2.8% is again estimated for the deep burnup element in 25-second measurement in the MPBR application 
under design. 

Figure [TOl further shows the profile of the measured net counting rate of ^^''Cs as a function of the input 
counting rate for CSl and CS2. It is shown that below 100 xlO^s"^ input rate, the net rate keeps constant, 
revealing a reliable extraction of the net counting rate of the EOI ^^^Cs . Again for CS2, corresponding to a 
deep burnup element in our design, the standard variance of the averaged net rate is within 3%. 



4 Summary 



In summary, we have studied experimentally the feasibility of utilizing an HPGe detector in the deter- 
mination of the fuel element burnup in a future bed-like reactor. In the multi source measurement, the 
passthrough curves are measured by varying the shaping and fiattop time of the spectrometer. The plateau 
of the passthrough rate covers the radiation background level from 100 to 200 xlO^s"^ with adequate energy 
resolution and allows to implement the activity measurement of the burnup indicator ^''^Cs within half minute. 
The energy resolution, reflected by the FWHM of the full energy peak of ^^''Cs in this measurement, allows to 
identify the peak of ^^^Cs with its various neighboring 7 rays despite a slight degradation of 20% at very high 
background level. The precision (Icr) of the net counting rate of the EOI ^^^Cs under various conditions is stud- 
ied. It is demonstrated that the subtraction of the background contributes significantly to the total precision of 
the net area of ^^^Cs peak. By fine tuning the background window and fixing the range of interest, the precision 
(Icr) is optimized to 2.8% with the typical incident radiation rate and ^^^Cs intensity that is relevant to the deep 
burnup situation in the future application. It is worth mentioning that because the background composition 
is much more complicated in the real reactor application than the experimental conditions presented here, a 
certain deterioration of the precision is expected and demands a detailed simulation (via MCNP or Geant4 for 
instance) of the real background irradiation. 
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Fig. 8. (Color online) Relative uncertainty of the net area as a function of Co irradiation rate for CSl 
and CS2. The curves denote the relative stand deviation of the net area distributions, and the symbols is 
calculated from Formulae ((3]) . 




Fig. 9. (Color online) Correlation between the relative uncertainty of the net rate and the net rate scaled 
with the irradiation rates 
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Fig. 10. (Color online) The net rate of the ^^^Cs measured at different ''"Co irradiation rate for CSl and CS2. 
The error bars denote the standard deviation of the net rate distribution in the repeated measurements. 
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